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A model to calculate nuclear responses considering short-range correlation effects is pre- 
sented. The model is applied to the study of electromagnetic responses induced by one-body 
operators. We calculate one- and two-nucleon emission responses and cross sections of the 
I ^^O and ^'^Ca nuclei in the quasi-elastic region, and we compare them with experimental 

data. 
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cn . 1 Introduction 

o 

' The study of atomic nuclei has been characterized in the last few years by the development 

. and the application of technologies able to deal with realistic nuclear interactions. In the 

nuclear structure jargon one calls realistic those interactions whose parameters have been 
fixed to reproduce the properties of two, and eventually three, nucleon systems. 

A large number of theories has been developed to tackle the problem of solving the 
O I Schrodinger equation with realistic interactions. Three or four nucleon systems are described 

with Faddeev [Q, ^ , Correlated Hyperspherical Harmonics Expansion Q ^] and Green Func- 
tion Monte Carlo techniques which solve the Schrodinger equation without approxima- 
tions. The Green Function Monte Carlo technique has also been applied, with great success, 
^ I to nuclei up to A=7 p. 

■ Mainly for computational reasons, the straightforward applications of the above theories 

to heavier nuclear systems is not affordable, therefore different approaches such as Cluster 
Monte Carlo 1^, Brueckner theory Exponential S |^, Correlated Basis Function (CBF) 
Q, have been developed with the aim of obtaining approximate, but still accurate, solutions 
of the many-body Schrodinger equation. 



Among these theories the CBF |10] has been well tested and widely applied in infinite 
nuclear systems like nuclear and neutron matter These studies have shown that, 

in the framework of the CBF theory, the Fermi Hypernetted Chain (FHNC) resummation 
technique with the Single Operator Chain (SOC) approximation provides solutions of the 
Schrodinger equation having an accuracy of about 1 MeV per nucleon on the binding energy. 



In this computational scheme the use of Argonne nucleon- nucleon potentials 15], together 



with Urbana three-body forces ||T^, provides equations of state whose minima are rather close 
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to the empirical one. In the same framework also the nuclear matter responses to external 
probes have been calculated obtaining satisfactory agreement with experimental inelastic 
electron scattering data |17, 15]. 



The success of the infinite nuclear matter results has lead one to apply the CBF theory to 
the description of the ground state of finite nuclear systems. Recently, ^^O and ^'^Ca ground 
states have been described within the FHNC-SOC computational scheme with interactions 
including tensor, spin-orbit terms and three-body forces |^] . 

In this paper we deal with the problem of describing excited states of finite nuclei within 
the CBF computational scheme. The model we present is inspired by the nuclear matter 
works of Refs. [l^. While in these works the CBF cluster expansion has been fully 
considered, in our model we retain only those terms containing a single correlation function. 
The validity of this truncation has been tested, for the nuclear matter charge response, in 
Ref. 11^. Because of the excellent agreement between the full calculation and our truncated 
model we felt confident enough to extend our model to include the current operators. 

Our model treats the short-range correlations but does not consider collective nuclear 
excitations. For this reason we have applied it to the description of the nuclear responses in 
the quasi-elastic region where the collective effects are neg ligible pOl, P]. 



In this paper we shall deal only with inclusive electron scattering data, and we consider 
excited states with both one and two particles in the continuum. The two nucleon emission 
is treated as a genuine short-range correlation effect, and we neglect the contribution of the 
two-body currents. 

After briefly reviewing in section § the approach of Ref. [0] , we present in section ^ our 
model, and we apply it in section |^ to inclusive electron scattering. The results obtained in 
the calculations of the ^^O and ^'^Ca quasi-elastic responses are presented in section ^ In 
section ^ we summarize our work and we draw our conclusions. 

2 Responses in a correlated theory 

The linear response of a many-body system to the perturbations induced by an external 
operator 0(q) is given by: 

5(q,cu) = --ImD(q,w), (1) 

TT 

with 

D{q, Lo) = (io I 0+(q) {H-EO-UJ + irj)-' 0(q) | ^o) , (2) 
where we have indicated with | the normalized eigenstates of the nuclear hamiltonian H: 

l^n> = --^- (3) 
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In the previous equations q and uj represent the momentum and energy transferred to the 
nucleus. 

Inserting in Eq. (^) a complete set of eigenvectors of H we obtain: 

n - ^0 - W + zr/ 

where we have defined: 

(^„|0(q)|^o) 
= 1 r • (5) 

Like in CBF theory we assume that the nuclear many-body ground state is described as 
a product of a correlation function G and a Slater determinant |$o) of a set of single particle 
wave functions occupying all, and only, the states below the Fermi surface: 

|^o) = G|$o). (6) 

Since the Slater determinant is already antisymmetrized with respect to the exchange of two 
nucleons, the correlation operator is given by a symmetrized product of two-body correlation 
operators: 



G(l,2...yl) = S 



Kj 



(7) 



where we have indicated with S the symmetrizer operator. 

In modern nuclear structure calculations with realistic microscopic interactions, the two- 
body correlation operator is taken as a sum of operator dependent correlation functions 



F^J = E r(-^.)Of,> (8) 

p=l,8 



where the involved operators are: 

Of~^'^ = [1, (Ti • O-j, Sij, (L • S)ij] (g) [1, Ti • Tj] , (9) 

with Sij = (3fjj • (Tj fjj • (Tj — (Ti ■ aj) indicating the tensor operator, and rij = |rj — rj| the 
distance between the positions of the particles i and j. 

The ground state wave function is obtained by minimizing the ground state energy with 
respect to variations of the correlation function and of the single particle basis. In the theory 
developed in Ref. the correlation operator G whose parameters have been fixed by the 
ground state minimization is used to generate the excited states of the system: 

\m^) = G\^n)- (10) 
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The mean-field excited states are obtained by making particle-hole excitations on |$o)- 
In order to use the cluster expansion techniques it is useful to rewrite the function ^n(q) 

as: 



^n(q) 



($„|G+0(q) G\<^ 



0/ 



($o|G+G|$o) 



{^o\G+G\^o) 

{<i>n\G+G\^n) 



(11) 



The two factors in Eq. (|llD are evaluated separately by expanding the numerator and the 
denominator in powers of the short-range function Fij — 1 (see Ref. [^] for a detailed 
presentation of the cluster expansion of Cn(q))- In both factors the denominators cancel the 
unlinked diagrams. 



3 The model 



Our model simplifies the calculation of the cluster expansion of Eq. (11) by retaining only 
those terms involving a single correlation line. This model has been inspired by the results of 
Refs. [27, 28| where the density and momentum distributions of doubly closed shell nuclei, 
calculated with a single correlation line model, have been shown to be rather similar to those 
obtained with complete CBF/FHNC calculations. 

In the present article we consider that the external operator 0(q) is a one-body operator 
and we use only scalar correlations. Considering Ofj = for p > 1 in Eq. (^), we can express 
the function ^„(q) as: 



?n(q) 



(^n|n.<, 40(q) U^<if^J\'^' 



i<j J^3 



(^n|0(q) Y{i<j{'^ + hij)\^ 
{'^o\Y\i<j{'^ + hij)\^o) 



('^oini<j(l + ^)l^o) 

(^n|^^<i(l + /i^J)l^n)_ 



(12) 



where we have considered that fij is a real scalar function, therefore commuting with the 
operator 0(q), and we have defined hij = ff^ — 1. 

The approximation of our model consists in retaining in Eq.([T^) only those diagrams 
having a single h function: 



^n(q) ^ ili^) = (^nl 0(q) E (1 + I^o)l 



(13) 



where the subindex L indicates that only the linked diagrams are considered. This result has 



been obtained using the same procedure used in Ref. |27] for the density distribution. The 
first step consists in making the full cluster expansion of numerators and denominators, and 
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Figure 1: Diagrams considered in the calculation of the one-particle one-hole responses. The 
dashed lines represent the correlation function h and the continuous oriented lines the single 
particle wave functions. The letter i,k,h indicates hole wave functions and p particle wave 
functions. A sum on the i and k indexes is understood. The black squares represent the 
point where the excitation operator is acting. 
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Figure 2: Diagrams considered in the calculation of the two-particle two- hole responses. The 
symbols have the same meaning as in Fig. |l]. 
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this allows the elimination of the unlinked diagrams. Only at this point do we truncate the 
obtained result by considering only the first order terms in hij. 

The terms contributing to C^(q) are presented as Mayer-like diagrams in Figs. 1 and 
2 for Ip-lh and 2p-2h excitations respectively. In each diagram, the black square indicates 
the coordinate where the excitation one-body operator 0(q) is acting, while the black dots 
indicate the other coordinates. The dashed line indicates the correlation function h, which 
operates on two-coordinates only, and the continuous oriented lines indicate the single particle 
wave functions (p. We have used the convention of considering entering into a point the wave 
functions of the state |$), while the wave functions of the state (^>| are exiting from the point. 
The letters k indicate single particle wave functions below the Fermi surface and imply a 
summation over them. With p and h we have labelled those particle and hole states whose 
quantum numbers characterize the full many-body excited state. 

Let us first consider the case when the nuclear final states have only one particle in the 
continuum. If we label with 1 the coordinate where the external operator 0(q) is acting, we 
can specify the in Eq. (13) as: 

A 

dpih(q) = (^ipihl 0(q) |$o) + (^ipihl 0(q) ^li I'^o) 

i>i 

A 

+ ($ipih|0(q) Y h^J\^o). (14) 

i<i<j 

The above expression shows that our model, in addition to the uncorrelated transitions 
represented by the the one-point diagram (1.1) in Fig. 1, generates also two- and three-point 
diagrams. The presence of both two- and three-point diagrams is necessary to have a correct 
normalization of the many body wave function. This can be seen by considering that, if 0(q) 
is the density operator, in the limit for g — > 0, when p = h, the sum on all the hole single 
particle wave functions should provide the proper number of nucleons A, as happens in Refs. 



1 27, 28 1 . In this limit, the diagrams containing a correlation line should not contribute. This 
can be seen in Fig. 1, by joining the p and h lines and considering that in the three-point 
diagrams the q ^ limit implies that the wave functions linked to the black square should 
be equal because of the orthonormality of the set of single particle wave functions. One can 
then observe that the contribution of the (2.1) and (2.3) diagrams is exactly canceled by that 
of the (3.2), (3.4) and (3.6) diagrams, and the contribution of the (2.2) and (2.4) diagrams 
is canceled by that of the (3.1), (3.3) and (3.5) diagrams. 

When the nuclear final state is characterized by two-particle two-hole excitations, the 
function ^^(q) is given by: 

A A 

^2p2h(q) = (^2p2h| 0(q) E ^li l^o) + (^2p2h|0(q) Y l^o) • (15) 

Kj l<i<j 
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As expected, the uncorrelated term does not appear, since a one-body operator cannot 
lead to a 2p-2h final state. The 4 two-point diagrams and the 12 three-point diagrams we 
consider are shown in Fig. 2. Also in this case the set of diagrams conserves the proper 
normalization. In analogy with the discussion done for the Ip-lh case one can see that, in 
the limit g — > and setting pi = hi and p2 = /12 the contribution of the diagrams of Fig. 2 
is zero, as expected. 



4 Electromagnetic excitations 



The response model presented in the previous sections has been applied to the description of 
inclusive electron scattering processes considering both one- and two- nucleon emission. 

In Plane Wave Born Approximation, the inclusive electron scattering cross section can 
be written as: 



dV 
dQ'de' 



q4 



+ tan^- 



2q2 



2p2h/ 



(16) 



where 9 is the scattering angle, q^ the four-momentum transfer satisfying the relation 
o;^ — q^, and o"m is the Mott cross section: 



CM 



acosie/2) 
2ei sin2(6l/2) 



(17) 



In the last equation a is the fine structure constant and ej is the initial energy of the electron. 

The electromagnetic responses Rl and Rt, depending on q = |q| and to, are obtained 
by evaluating the responses S{q,u!) in Eq. (Q) for charge and current operators respectively, 
and multiplying them by the electromagnetic nucleon form factors. In Eq. (16) the Ip-lh 
and 2p-2h responses do not interfere with each other because they produce two different final 
states. 

We have calculated the longitudinal responses using the charge operator: 



0(q) ^ p(q) 



dV e''!-'' pir) 



with 



= 2^ ^(^-^k) 



(18) 



(19) 



k=l 
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where T3(/c) is the third component of the isospin of the k-th nucleon. 

The transverse responses have been calculated considering convection and magnetization 
currents: 

0(q) ^ J(q) = j'^(q) +j-(q) = J dV e'^ '^ [j^(r) +j"^(r)], (20) 
with ^ 

= E ^ [^(r _ r,) Vfc + Vfc 6{r - r,)] (21) 

and 

j"^(r) = E {^^'-^Y^ + ^^^^?^) X - '^(^) • (22) 

In the above equation we have indicated with Mfc and cr(k) the mass and Pauli spin matri- 
ces corresponding to the A:-th nucleon, and with fip (^at) the proton (neutron) anomalous 
magnetic moment. 

Since in the quasi-elastic peak the contribution of j"^ is small compared to that of (see 
the discussion of Fig. 5 in the next section) we have considered the convection current only 
in the uncorrelated one-body term (1.1) of Fig. 1. 

We restrict our study to the investigation of doubly closed shell nuclei and we suppose that 
the target nucleus makes a transition from its ground state to an excited state characterized 
by the total angular momentum J, its projection M on the quantization axis, and the parity 
n = ±1. The responses can be written as: 

oo 

RLiq,u;) = E \{^%\MjMiq)\^lo)\'SiEj-Eo-u;) (23) 
J,M 

and 

oo 

RTiq,u;) = E {K^Jm|T|m(9)I^Jo)I' + \{^%\Tjiliq)\^lo)\'} S{Ej - Eo - u;) , (24) 
J,M 

where we have not indicated the parity of the ground state which is always positive in our 
case. The Mjm and Tjm operators used in the previous equations are obtained by making 
a multipole expansion of Eqs. (|^) and (pO[). For the charge we have: 

Mjm (q) = Jd'rjjiqr) Yjm (r) p(r) , (25) 

where we have indicated with f the 9 and (j) angles characterizing the vector r in polar 
coordinates, with Yjm the spherical harmonics and with jj the spherical Bessel functions. 
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For the current we should distinguish between the electric excitations (E) with natural 
parity 11 = ( — 1)"^ and the magnetic excitations (M) with unnatural parity 11 = (— l)"^"^^: 

TfM = - fd\{Vx [jj{qr) Y)S(f-)] } • J(r) (26) 



and 



T^Miq) = J dVij(gr) Yf^(f) • J(r) . (27) 



where we have used the symbol Yjj to indicate the vector spherical harmonics p 



We calculate the responses using the expressions (||), (^) and (|l^), and we rewrite Eqs. 
(H) and (|§ as: 

oo 

RLiq,uj) = i7rJ2\{<^'J\\Mjiq)J2i^ + hij)\\'^h)\''SiEj-Eo-uj) (28) 
J i<j 

and 

oo I 

RT{q,u;) = 4vr^ KcI>n||rJ(g)^(l + /i,,)ll<J>J)P 
J [ i<j 

+ Tfuiq) E (1 + \ KEj -Eo-u). (29) 

i<j J 

The above expressions have been obtained by applying the Wigner-Eckart theorem, there- 
fore the symbol " || " indicates that the angular part should be calculated by considering only 
the reduced matrix elements. 

The calculations of the transition matrix elements are carried out by performing a mul- 
tipole expansion of the correlation function h: 

oo 

hij = h{rij) = h{r i,rj, cos 8ij) = E hL{ri,rj) PL{cos9ij) , (30) 

L=o 

where Pl are Legendre polynomials. Because of the completeness and the orthogonality of 
the Legendre polynomials, the hLiri,rj) terms can be evaluated as: 

2L + 1 

hL{ri,rj) = — ^ — J ^d{cos 9ij) Pl{cos 9ij)h{ri,rj, cos 9ij) . (31) 



Eqs. ( p8[ ) and ( p9[ ) show that the final state quantum numbers are determined by the 
quantum numbers of the uncorrelated many body state l^j^v/) which is built as a Slater 
determinant of single particle wave functions of the form: 

Mr) = Rii,{r) E i^f^Ujm) Yi^i^) xs Xt , (32) 
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where Xs and xt are the spin and isospin wave functions respectively and {Ifi^sljm) a Clebsch- 
Gordan coefficient. 

We treat now separately the one and two nucleon emission cases. In the first one, once the 
hole is fixed, energy conservation defines the value of the emitted particle energy: gp = uj + eh- 
The angular momentum coupled Slater determinant is related to the uncoupled one by the 
relation: 

\^%{lvlh))= {-l)"{jpm,jf,-m,\JM)\^l,^), (33) 

mpirih 

where jp and jh are the particle and hole angular momenta and nip and m/j their third 
components. We insert the expression (53) into Eq. (14) and apply the Wigner-Eckart 



theorem. Using the orthornormality of the single particle wave functions we obtain the 
following expression: 

A A 
^ipM = iP\\ Oj(q) \\h) + E iPj\\ Oj(q) hu \\hj) + E {pjk\\ Oj(q) h^sWhjk) . (34) 

j j,k 

In this equation we have labelled with 1 the coordinate where the electromagnetic one-body 
operator is acting. The second and the third terms of Eq. (^) produce the two- and three- 
point diagrams of Fig. 1 respectively. 

The multipole expanded expression of a generic one-body operator is formed by the 
product of a term depending on q and ri multiplied by a term depending only on the angular 
coordinates: 

Ojm{^) = Fj{qn)OjM{h). (35) 

In the case we are interested this is obtained by inserting the expressions (|l9|), ( |2lD and 
(H) into Eqs. @) and (H). 

We present in the following the expressions of the terms of Eq. ( p^ ) using the notation 
of Ref. [^] to indicate the 3j, 6j and 9j symbols of Wigner. For the one-point diagram we 
obtain: 

T^i^P^^) = o;/^ J dn rl R;{n) Fj{qn) Rhin) , (36) 

where we have defined: 

={lalja\\OK{h)\\hlj/3)- (37) 

The expressions we have obtained for the two-point diagrams are: 

dJp^^) = i7ru;'^j,ldnrlR;in)Fjiqn)Rf,in)ni'\n), (38) 

Z)Sr^ = 4vr(-l)^W^5:5,,J-l)^^7i{ f f ■] (39) 
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J dnrlR;{n)Fj{qn)R,{n)nf\n)^ 



J dri rl R*{n) Fj{qn)Ri{n) 'Hf\n) 



Dfr) = 4vr(-l)>+^'^ V<^,. (-l)^^X^J f ^ (41) 



iLK ^ I 



Jh Jp Ji 



In the previous equations, we have used the notation a = \/2a + 1 and the following definition 

of lafi- 

1 1 

7a/3 = (^a-ja II i^A(?) 1^/3 2^/3) 

= -l=(-iy/^+A-le(/„ + /^ + A)jJ^A ^j, (42) 

with ^(a)=l if a is even and ^(a)=0 otherwise. The other symbols used are: 

^I'p"" = («ll [ydOi. <Ai) ® Oj{l)f 11/3) (43) 

and the integrals 

n^f\ri)= j dr2rlRl{r2)hL{n,r2)Rp{r2) (44) 

and, in Eq. (|3|), 

^f' = i E ^0 (n) = / dr2 ri h,{n,r,) p{r,) , (45) 

where 

is the nuclear density. 

For the three-point diagrams we have obtained the following expressions: 

DiT^ = Y.^t,,tju,t,s,^,^aip + ik + L)aip + k)l ( \' \ \\ m 

ikL V 2 2 ^ / 
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^S''^ = 47r 5: 5t.t, 5,.,, + y J dri rf i?*(ri) Fj{qn) Rn{n)J^'^-^^''^ , (48) 



2 



ikh 

u^i, I dnrlR;{n)Fj{qn)Rkin)j!;''^-'^'''^ 
= 47rY,St,,tJ,,,,ak + lh)u^',i JdnrlR;{n)Fj{qn)R,{n)jt^^^^^^ (50) 



^.i / dnrlR*in)Fjiqn)R,in)jt^'^''^ 

^ir'^ = 4vr<5,^,,7/.iE^*.,*.7/.a;,i/dnr?i?|(r0F,(,r0j}^^]^['^^^ (52) 



ik 

In these equations we have introduced the symbols: 



Jt'^'^''^ = J dri rl Kin) R,{n) ^[^^'(n) (53) 

^ Ur^rlRl{n)Rp{n)nf{n). (54) 



and 



The equations presented above have been obtained for a generic one-body transition 
operator of the form of Eq. (^). The calculation of the electromagnetic response functions 
continues by inserting in the above equations the explicit expressions for the charge and 
current operators, Eqs. (|25|)-(p7[). The final expressions are given in Appendix A. 

In the case of two nucleon emission the total angular momentum and parity of the nuclear 
final state are given by the combination of the four angular momenta and the parities of the 
particle and hole states involved in the excitation. We have chosen to couple the angular 
momenta of the two hole states and of the two particle states first, and then, with another 
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recoupling, to obtain the total angular momentum of the 2p-2h excited state: 

|$V(2p2h)) = J2 {JpMpJh-Mh\JM) J2 {jp,mpJp^mp^\JpMp) 

E (-l)^-'>i+^-'^2--'H— - mnjh, - mh,\JhMH) |1>?p2h) • (55) 

Inserting Eq. (p5|) into Eq. (15) and applying the Wigner-Eckart theorem we obtain the 
following expression: 

A 

^2p2h(q) = {PmW Oj(q) hu \\h1h2) + E {PiP2k\\ Oj(q) /i23 \h1h2k) , (56) 



where the coordinate labels have the same meaning as in Eq. (p^). In the 2p-2h excitation 
the energy conservation does not uniquely define the energies of the emitted particles since 

w = + - e/ii - et^ , (57) 

for this reason in Eq. ( [55| ) an integration on the energy of the particle p2 is implied. 

It is evident from Fig. 2 that the various diagrams can be related one to the other by 
changing the labels of the hole and of the particle states involved. We use the following 
symmetry relations to calculate the contribution of the various diagrams to the response: 

D^2f''\pihiP2h2] = {-1^^+^''^+'- D^^f''\p2hmh2], (58) 

D!^f''\pihiP2h2] = (-l)M+iP2+^P (-l)^'^i+J'^2+^hDgP2h)[p2/i2Pi/li]. 

The explicit expression of D^f'^^^ matrix element for a generic one-body operator is given 
by: 

^ I jpi jhi K 



D^2.f''\pihip2h2] = 47r5t^^,t,^ JpJ,E(-l)''T^7p2/.2 jh, L \ (59) 

LK ^ [ Jp Jh J 

I dn rl r;^ {n)u;l;f^f R,, (n) n^''^'^ (n) . 

The symmetry relations we used to evaluate the the three-point diagrams are: 
D^f"^ \p1h1p2h2] = {-iy'n+^'^2+J'^ D^^f"^ [pih2P2hi] , 

Dif''^\pihiP2h2] = {-iy^-+'^'^+'^Dif''\p2hip^h2], 
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Dif''\p,h,p2h2] 
D^^f''^\pihiP2h2] 



D 



if^^ \p1h1P2h2] 



(_l)ipi+iP2+-^P ^_iyh,+h^+JH D^^f^\p2h2Pihi] 
{-iy^i+^'^2+J^D^^f''\pih2P2h,], 

{-iyvr+J.2+J,D^^f^^\p2hiPih2], 

(_iyPi+iP2+-^p ^_iyh,+jH^+JH D^^f^\p2h2Pihi] 
(_l)ipi+iP2+-^P ^_iyH,+jH^+JH D^^f^\p2h2Pihi] 
(_l)ipi+iP2+-^p (-iyh,+h2+Jh Di^f^^\p2h2Pihi] 



Dif2^\pihip2h2] = 
The corresponding matrix elements for a generic one-body operator are: 



D^^f''\p,h,p2h2] = ^ndt^^,,^ J, Y: 5u,t,, 



iL 



L2 



^P2h2 ^ihi 



jp2 '^h ji 
jh^ L jh2 



(60) 



(61) 



Jp 'Jh J 

Ji Jpi Jp2 



Dff''\ihiP2h2] 



^^K^t,^ {-ly^^^^'-^^'^'jpJh E K,t,, ^7p,h,lp2^ 

iL ^ 



(62) 



ji jhi Jp 
jpi jp2 ^ 



Jp 'Jh J 

Jh2 Ji Jhi 



[P\h-i];\p2i] 



Dif'\p,h,p2h2] 



^^K2^t,^K2,j>n i-iy^^^^^-^'^-^'-- 



r \ (63) 

Jh2 Jpi Jp2 I 
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Also in this case the contribution of the electromagnetic charge and current operator to 



Jp Jh J 

Jh2 Jpi Jp2 



the two nucleon emission responses has to be calculated by inserting Eqs. (|25D-(P7|) in the 
above equations. The explicit expressions we have obtained are given in Appendix B. 
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5 Specific applications 



The model described above has been apphed to evaluate electromagnetic responses and cross 
sections in the quasi-elastic region for the doubly magic nuclei ^^O and ^'^Ca. 

The two-body correlation function and the single particle basis are the two inputs nec- 
essary to build the many-body wavefunctions used in the calculation of the responses. In 
our calculations the set of single particle wave functions have been generated by a spherical 
Woods-Saxon potential of the form: 

V{t) = VoF{r) + f — y ^ ^F{r) \ ct + Vc{r) , (65) 

where m.,^ is the pion mass, Vc{r) the Coulomb potential produced by a homogeneously 
charged sphere of radius R, and 

F(r) = ^ ■ (66) 

The same parameterization of the potential has been used for both bound and continuum 
waves. 

As we have already anticipated, the other input of our approach, the two-body correlation 
function, is a purely scalar function. Single particle waves and correlation functions are not 
independent quantities. In the theoretical framework of the CBF theory they are fixed by the 
minimization of the ground state energy functional. In our calculations we use wave functions 
and correlations taken from CBF-FHNC calculations of ^^O and ^'^Ca ground states done with 
different hamiltonians. 

A first set of many-body wavefunctions has obtained in Ref. [23| with the semi-realistic 
S3 interaction of Afnan and Tang. Specifically, we used the Woods-Saxon parameters of Tab. 
5 of that reference together with the correlations labelled Euler in Tab. 6. We should recall 
that in that calculation the Woods-Saxon parameterizations was taken from the literature 
and the energy minimum was obtained by changing only the correlation function. 

The second set of wave functions has been taken from Ref.|^] where a hamiltonian con- 
taining a V8' Argonne potential, which is a truncated version of the V18 interaction, together 
with the Urbana IX three-body interaction has been used. Our calculations have been done 
with the wavefunctions of Tab. V of that reference. Also in this case the Woods-Saxon 
parameters have been kept fixed and the energy minimum was found by changing only the 
correlation function. We should remark however that the energy minima found in these re- 
stricted calculations are only 7% larger than the those found by the full minimization. In 
our calculations we used only the scalar terms of the correlation functions of Ref. which 
consider the channels of Eq. (|^) up to p = 6. 

Henceforth we shall label S3 and V8 all the quantities related to the wavefunctions of 
Ref.[p^] and P] respectively. We show in Fig. 3 the correlation functions /(r) and we give in 
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Table 1: Parameters of the Woods-Saxon potentials used to generate the set of single 
particle wave functions. 
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0.53 


0.5 


0.53 



Table 2: Single particle energies, in MeV, of the O nucleus obtained with the two 
chosen potentials. We have indicated with tt and v proton and neutron states respec- 
tively. 
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Table 3: Single particle energies, in MeV, of the ^°Ca nucleus obtained with the two 
chosen potentials. We have indicated with tt and v proton and neutron states respec- 
tively. 







4"Ca 






V8 


S3 


c\xp 


TT lsl/2 


-32.5 


-36.0 




lp3/2 


-24.0 


-26.5 




lpl/2 


-24.0 


-23.0 




ld5/2 


-14.0 


-16.2 


-14.7 


ld3/2 


-14.0 


-8.7 


-8.3 


2sl/2 


-11.0 


-10.3 


-10.3 


V lsl/2 


-40.4 


-42.6 




lp3/2 


-31.5 


-32.4 




lpl/2 


-31.5 


-29.8 




ld5/2 


-21.1 


-21.4 


-18.2 


ld3/2 


-21.1 


-15.7 


-15.6 


2sl/2 


-18.2 


-16.4 


-21.6 



18 



1.20 
1.00 



0.80 - 



0.40 
0.20 
0.00 
1.20 
1.00 



1 1 1 1 1 


yv 
// 
// 
// 




if 


S3 




V8 


// 
// 
// 
/ / 




y / 


16^ 


1 1 


1 1 1 



0.0 0.5 



0.80 - 



0.40 



0.00 



1.0 



1.5 2.0 2.5 3.0 



1 1 1 1 1 


_ A 


S3 




// 
// 
/ / 


V8 




/ / 
/ / 
✓ / 


40^ 
La 




1 1 


1 1 1 





0.0 0.5 1.0 1.5 2.0 2.5 3.0 
r[fm] 

Figure 3: Correlation functions used in the calculations. The label S3 refer to the FHNC 



calculations of Ref. |23] while V8 to those of Ref. 
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Tab. |l] the values of the parameters of the Woods-Saxon well. To complete the information 
on the input, we present in Tab. § and |^ the values of the single particle energies and we 
compare in Fig. 4 the MF and the correlated charge density distributions with the empirical 



densities of Ref. |24]. The charge densities have been obtained by folding the pointlike proton 



densities with the electromagnetic nucleon form factor of Ref. |25|. The same nucleon form 
factor has been used in the calculation of the responses. 

Since the convection current contribution is already small at the MF level, as we show in 
Fig. 5, we did not consider the correlated terms produced by this current. In any case the 
convection current is included in the mean-field (MF) term. 

5.1 One- nucleon emission 

In Figs. 6, 7, 8, and 9 we present the electromagnetic Ip-lh responses calculated for the two 
nuclei considered for three different values of the the momentum transfer. 

In these figures the full lines indicate the MF results. The dotted lines have been obtained 
by adding to the MF responses the contribution of two-point diagrams of Fig. 1, while the 
dashed lines show the results of the calculations when all the diagrams have been considered. 

It is evident that the short-range correlations produce small effects on the inclusive re- 
sponses. In the peak positions we found maximum relative variations with respect to the 
MF responses of 1.7% in the longitudinal and of 2.2% in the transverse responses. These 
values are within the range of the uncertainty produced by the different choices of nucleon 
form factors ||2^ . In spite of this we made a detailed investigation of the correlation effects 
because they can be relevant in the study of exclusive experiments. 

In order to have a better view of the correlation effects we show in Figs. 10 and 11 the 
difference between the correlated and MF responses: 

Ai?L,T('7, oj) = RL,T{q, io) - Rf^iq, ^) • (67) 

In these figures we have used the convention of indicating with the dotted lines the results 
obtained with the two-point diagrams only, and with the dashed line those obtained with all 
the diagrams. 

The inclusion of the three-point diagrams reduces the effects produced by the two-point 
diagrams alone. This fact is more evident in the longitudinal responses and it is consistent 



with the findings of Refs. |27, 28] where ground state charge and momentum distributions 
have been studied. As discussed in sect. ^, the proper normalization of the many-body wave 
functions is obtained because in the limit for q ^ the three-point diagrams cancel exactly 
the contribution of the two-point diagrams. This cancellation is present also for large values 
of the momentum transfer, but it is not exact any more. It is interesting to notice that the 
final correlated longitudinal responses are always smaller than the MF responses with the 
only exception of the S3 result in ^'^Ca at 300 MeV/c. 
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Figure 5: MF transverse responses calculated with the S3 Woods-Saxon potentials consid- 
ering magnetization current only (dashed lines) and magnetization plus convection current 
(full lines). The numbers in the figure indicates the values of the momentum transfer in 
MeV/c units. 
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Figure 6: Longitudinal and transverse Ip-lh responses in ^^O for three different values of 
the momentum transfer. The calculations have been done with the V8 wave functions and 
correlations. The full lines show the MF responses, the dotted lines have been obtained with 
the inclusion of the two-point diagrams of Fig. |l|, while the dashed lines have been obtained 
considering all the diagrams of that figure. 
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Figure 7: The same as in Fig. with the S3 input. 
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Figure 9: The same as in Fig. |^ for ^°Ca. 
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Figure 10: Differences between correlated and MF responses as defined in Eq. (|67|) for ^^O. 
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Figure 11: The same as Fig. 11^ for "^^Ca. 
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The behavior of the correlations on the transverse responses is rather different: the cor- 
related responses are always larger than the MF ones. The reason of this behavior has to be 
ascribed to the diagram (2.3) of Fig. 1 which does not contributes in the transverse response 
as indicated in Appendix A. In this diagram the exchanged virtual photon strikes a nucleon 
which remains bound but through the correlation it excites a p-h pair producing the nucleon 
emission. Since in our calculations the correlation is a scalar function, only longitudinally 
polarized virtual photons can induce this type of excitation. 

To study the effect of this diagram we set it artificially to zero also in the longitudinal 
channel. The comparison with the correctly calculated longitudinal response is shown in Fig. 
12 for the ^^O response at q = 400 MeV/c calculated with the S3 wavefunctions. In the 
upper panel we show the responses, and in the lower one the ARi of Eq. (|67|). The (2.3) 
diagram lowers the response and produces the discussed difference between longitudinal and 
transverse responses. The curves without the contribution of the (2.3) diagram in the lower 
panel are very similar to the analogous ones of Fig. 10 for the transverse responses. 

Another interesting feature observed in Figs. 10 and 11 is the different behavior shown 
by the S3 and V8 correlations. The two-point contributions of the S3 correlation are always 
positive, and the S3 total contributions are always larger than those of the V8. 

The different features characterizing the two correlation functions are shown in Fig. 3. 
A first difference concerns the behavior at r = where the S3 correlations assume larger 
values than the V8 ones. The second difference is related to the behavior at r values around 
1-1.5 fm, where the S3 correlation functions show a sizable overshooting with respect to the 
asymptotic value of 1. In order to study the effects of these two features on the response 
functions, we performed calculations using correlation functions of the form: 

/(r) = 1 -Aexpl-Br"^] + Cexp[-D{r - rof] . (68) 

This expression allowed us to disentangle and magnify the characteristics of the correlation 
function under investigation by changing the values of the parameters. The values used in our 
calculations are given in Table ^ The correlations gl, g2 and g3, shown in the upper panel 
of Fig. 13, are composed only by a single gaussian and have been constructed to change only 
the behavior at r = 0. The other correlations, shown in the lower panel of Fig. 13, have the 
same behavior at r = but have different overshooting properties at intermediate distances. 
With these correlations we calculated the ^^O responses at q=400 MeV/c using the S3 set of 
single particle wave functions. The results of these calculations are summarized in Fig. 14 
where the ARl,t defined in Eq. (|67|) are shown as functions of the excitation energy. 

The two upper panels of the figure show the effect produced by changing the correlation 
depth at r = 0. The global effect is a reduction of the MF responses, more evident in the 
longitudinal one. The reduction is larger for deeper correlations. The other two panels of the 
figure show the effects of the overshooting. There is a slight depletion of the longitudinal MF 
responses and an enhancement of the transverse ones. The effects are magnified by shifting 
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Figure 12: Upper panel: ^^O longitudinal responses calculated with the S3 input. The curves 
indicated with NO (2.3) have been obtained without the (2.3) diagram of Fig. |. The other 
curves correspond to those shown in Fig. |^. Lower panel: difference between the correlated 
and MF responses. Also in this figure the dotted lines indicate the results obtained by using 
only the two-point diagrams, while the dashed lines show the results of the full calculation. 
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Figure 13: Correlation obtained with the expression ( pq ) for the various parameterizations 
given in Tab. 0. 
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Table 4: Parameters used for the correlation given in Eq. (|6^). 



gl g2 g3 g4 g5 g6 

A 0.8 0.6 0.4 0.8 0.8 0.8 

B [fm-2] 2.2 2.0 1.9 2.2 2.2 2.2 

C 0.0 0.0 0.0 0.04 0.04 0.04 

D [fm-2] 0.0 0.0 0.0 3.0 3.0 3.0 

ro [fm] 0.0 0.0 0.0 1.5 2.0 2.5 



the overshooting at larger distances. This figure clearly shows the role of the three-point 
diagrams which always reduce the effects of the correlations. 

To compare of our results with experimental data it is necessary to consider the effects of 
the rescattering of the emitted nucleon with the residual nucleus . We describe these Final 
State Interactions (FSI) effects with the method developed in Refs. |20, The method 
shows that, under certain approximations, it is possible to consider the FSI by folding the 
bare responses with a smoothing function: 



R''^\q,io) = — dER{q,E) 



M* 



M 



p E, —rUJ ] + p[E, — —UJ 



M* 



M 



(69) 



where we have indicated with M and M* the free and effective nucleon masses respectively, 
and where the function p{E, to) describing the effects of FSI is: 



p{E,oj) 



1 



(70) 



27r[E-iu- A{lu)Y + [r{uj)/2] 
The A{uj) and r(a;) functions of the above equation are connected by a dispersion relation: 

1 



2tt 



-V 



duo — ■ 



UJ' 



■ UJ 



(71) 



where V indicates the principal value of the integral. The function T{uj) is related to the 
imaginary part of the single particle self-energy. We have obtained r(a;) as average of the 
single particle energy width 7(1-1;): 



1 /""^ 

T{uj) = - de[-f{e + uj) + -i{e - uj)] . 
w Jo 



(72) 



The single particle energy has not been calculated with a microscopic model, but, following 
Ref.[^], has been parametrized as follows to reproduce empirical values: 



7(e) 



+ 62 



%), 



(73) 
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with a = 11.5 and b = 18. 

The effective nucleon mass takes into account non-locahty effects of the mean-field which 
cannot be neglected. We determined the value of the effective nucleon mass by using the 
expressions proposed in the polarization potential model of Ref. 

Ml = I (74) 

with 

^2 ^ 

M/M, 

AC/=50 MeV and = 0.8M. 

In Fig. 15 we compare our ^'^Ca longitudinal and transverse responses with the data 



2 _ 2MAU 
^° " M/M* - 1 ' 



of Refs. 1 32] and |33]. Dashed (full) lines show the results of the S3 (V8) calculations. 
Both results are in a reasonable agreement with the data. The differences between the two 
calculations are mainly produced by the differences in the mean-field basis rather than in the 
correlation functions. 

For the nucleus ^^O, Rosenbluth separated data for longitudinal and transverse responses 
are not available. In Fig. 16 we compare our results with the cross section data taken at 
Frascati These data have been measured for fixed values of the incident electron energy 
and of the scattering angle and, therefore, the value of the momentum transfer changes at 
every excitation energy. Since the data have been taken at the scattering angle of 32°, the 
average momentum value for the cross section at e=700 MeV is about 400 MeV/c and that 
at e=1200 MeV is about 650 MeV/c. 

The folding method used to consider FSI is thought for fixed values of the momentum 
transfer. Every cross section point has been calculated by fixing the value of q, calculating 
the response on a wide energy range (up to 400 MeV), folding the response with Ed. (|6g|), 
and extracting from this output the single result at the appropriate excitation energy value. 
The responses calculated in this way have been inserted in Eq. (p!6|). 

Also in Fig. 16 the full (dashed) lines show the results obtained with the V8 (S3) input. 
The agreement with the data is satisfactory, especially for the cross sections measured at 
higher energies. The differences between the S3 and V8 results are even smaller than in the 
case of ^'^Ca and the data cannot discriminate between them. 



5.2 Two-nucleon emission 

The study of the two-nucleon emission responses is limited by the heavy computational load. 
For this reason we have restricted our investigation to the ^^O nucleus only. 

For these calculations we used the same configuration space considered in the one-nucleon 
emission calculations. This includes all the protons and neutrons waves with I < 14. Within 
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this configuration space all the excitation multipoles compatible with angular momentum 
coupling of the 2p-2h excitation pairs have been considered (see Eq. (|55|)). 

We found that the contribution of the various partial waves becomes smaller the higher 
is the angular momentum value. Specifically, we have seen that for q ~ 600 MeV/c the 
contribution of partial waves with / > 4 is about the 2% of the total value at an excitation 
energy of 100 MeV and becomes about the 40% at the excitation energy of 300 MeV. These 
percentiles remain roughly constant in the momentum transfer range between 400 and 600 
MeV/c. In order to reduce the computational time we performed a complete calculation for 
the partial waves with / < 4 while we made a statistical evaluation of the contribution of 
the higher partial waves. This statistical evaluation consist in calculating a randomly chosen 
10% of the 2p-2h matrix elements, and in rescaling the obtained results proportionally to 
the total number of the matrix elements to be calculated. This simplification reduced by 
2/3 the computational time. We have tested the reliability of this simplifying technique by 
calculating all the 2p-2h matrix elements at the excitation energies of 100, 200 and 300 MeV. 
In all the cases investigated we found differences with respect to the simplified calculations 
smaller than 1%. 

The calculations of the two-nucleon emission responses have been done for both correla- 
tion functions at excitation energies of 100, 140, 200, 240 and 300 MeV for the kinematical 
conditions used to evaluate the inclusive cross sections of Fig. 16 at e = 700 and 1080 MeV. 
The values of the responses obtained in these calculations are shown in Fig. 17 by the squares 
(S3) and by the circles (V8). 

A direct comparison with the results of Ref. [35| is not straightforward, since the calcula- 
tions of that reference have been done on the ^^C nucleus, and obviously, with different inputs. 
Considering these differences we can say that the order of magnitude of the two calculations 
is the same. In the two calculations it is however remarkably different the contribution of 
the three-point diagrams. We found a maximum contribution to Rl of about 8% and even 
smaller, 1%, for the transverse response. In the longitudinal response the interference be- 
tween two- and three-point diagrams is always destructive, therefore the final response is 
always smaller than that of the two-point diagrams alone. The contradiction with the find- 
ings of Ref. 1^], where the three-point diagrams were giving an important contribution to 
the response, was produced by a computer error we found in our previous calculations. 

The results of Fig. 17 show a large sensitivity to the choice of the correlation. This 
is not surprising since in these calculations the two-nucleon emission is produced only by 
means of the short range correlation. We have also verified the sensitivity to the choice of the 
mean- field by calculating responses with V8 correlation and S3 mean- field for e=1080 MeV. 
The results obtained in this way are shown in Fig. 17 by the triangles. The changes of the 
mean- field wave functions affect considerably the final result. 

The two-nucleon emission responses have been used to calculate the inclusive cross section. 
The total responses have been obtained by adding one- and two-nucleon emission responses. 
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To make this sum we have interpolated the results of Fig. 17 in order to obtain the values 
of the 2p-2h responses at the desired excitation energies. We have used linear and quadratic 
interpolations and their results are shown in Fig. 17 by the dashed and full lines respectively. 
These interpolations as well as more sophisticated interpolation techniques we have also used 
produce results which are the same within the accuracy of the calculations. For the evaluation 
of the cross section we used the results of the quadratic interpolation. 

The comparison of the results of Fig. 17 with those of Figs. 6 and 7 shows that the 
2p-2h responses at the peak energy are two orders of magnitude smaller than the Ip-lh ones. 
The contribution of two-nucleon emission becomes relevant only at higher energies, in the 
dip region. As example, we show in Fig. 18 the contribution of the various V8 responses in 
the dip region for the two kinematical cases we have considered. The dashed lines represent 
the tail of the Ip-lh responses while the dashed-dotted lines the contribution of the 2p-2h 
responses. The full lines are the sum of the two. It is interesting to notice the different 
behavior of the responses for the cases at 700 and 1080 MeV. In the first case, for excitation 
energies above 220 MeV the 2p-2h responses become larger than the Ip-lh ones. For e=1080 
MeV the Ip-lh responses remain always larger than the 2p-2h ones. This different behavior 
is produced by the larger value of the momentum transverse in this last case which shifts the 
Ip-lh responses to higher energies. For this reason the tail of this response has not yet died 
out as in the 700 MeV case. 

The calculation of the cross section requires the inclusion of the FSI effects. We have 
treated them with the folding model previously described. There arc ambiguities in our 
procedure of including the FSI interaction. We could apply the folding procedure to the 
total responses or separately to the Ip-lh and 2p-2h responses and then sum the two folded 
responses. The two procedures gave results differing only by few percent. We show in Fig. 19 
the tail of the calculated cross sections. The dashed lines present the pure Ip-lh results while 
the full lines have been obtained by including also the 2p-2h contribution. The effects of the 
2p-2h responses which appeared in Fig. 18 are washed out by the inclusion of the FSI, whose 
effects are much larger than those produced by the short range correlations. The inclusion of 
two-nucleon emission responses does not modify the comparison with the experimental data 
shown in Fig. 16. 

6 Conclusions 

In this paper we have presented a model to evaluate the effects of the short-range correlations 
in the excitation of the one-body nuclear responses. Our model takes into account both two- 
and three-body diagrams where the short-range correlations appear at the first order. We 
have shown that both kinds of diagrams are necessary to have properly normalized wave 
functions. 
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Figure 18: Longitudinal and transverse responses of in the dip region calculated with 
V8 input. The dashed lines show the Ip-lh responses, the dashed dotted lines the 2p-2h 
responses and the full lines show the total responses. 
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Figure 19: Inclusive cross sections on ^^O for two values of the electron incoming energies 
and a scattering angle of 32°. The full lines have been obtained including both Ip-lh and 
2p-2h contribution while the dashed lines show the Ip-lh one. 
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As input of our calculations we have used scalar (Jastrow) correlation functions and single 
particle wave functions fixed in Refs. |23| and Q by a minimization of the nuclear ground 
state energy. 

The model has been applied to evaluate the inclusive electromagnetic longitudinal and 
transverse responses in the quasi-elastic region for the ^^O and ^'^Ca nuclei. We have con- 
sidered both one- and two-nucleon emission channels. The longitudinal response has been 
calculated using the charge operator. In the transverse response we have considered the 
magnetization current and we have included the convection current only in the uncorrelated 
diagram. 

Our results show that the effect of the correlations on the one-particle emission responses 
is small, certainly within the range of uncertainty related to the different choices of the 
electromagnetic nucleon form factors. In spite of this, we made a detailed investigation 
of short-range correlations effects in our model because they can be relevant in exclusive 
processes. 

The three-point diagrams, usually neglected in the literature, interfere destructively with 
the two-point diagrams. For this reason any effect produced by the two-point diagrams is 
always reduced when the three-point ones are included. In certain situations this reduction 
is of the order of 80%. 

Our calculations show that the correlations act differently on longitudinal and transverse 
responses. Specifically, we found that the longitudinal responses are lowered while the trans- 
verse ones are enhanced. This happens because some of the three-point diagrams do not 
contribute in the longitudinal response. 

We have calculated the inclusive cross sections in ^^O to compare our results with the 
data of Ref. |34| . In ^''Ca the comparison has been done directly with the empirical responses 
of Ref. [^ ]. In both cases the agreement with the data is quite satisfactory independently 
of the inclusion of the short-range correlations. 

The two-nucleon emission responses have been calculated only for the ^^O nucleus. Be- 
cause of the heavy computational load we have simplified the calculation by evaluating the 
most important matrix elements and providing a statistical estimation of the others. 

As in the one-nucleon emission case, also in the two-nucleon emission responses, the two- 
and three-point diagrams produce effects with opposite sign. However, in this case we found 
that the contribution of the three-point diagrams is at most of the order of a few percent of 
that of the two-point diagrams. 

In the quasi-elastic region, our calculations show that the 2p-2h responses are two order 
of magnitude smaller than the Ip-lh ones. 

Summarizing our results, we should say that a proper normalization of the many body 
wave function in a first order expansion model for the correlations, requires the inclusion of 
the three-body terms. The contribution of these terms always reduces the effect produced by 
the two-point diagrams. This effect is not negligible in the one-nucleon emission responses, 
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but it becomes rather small in the two-nucleon emission. 

The use of different correlations, together with different set of single particle wave func- 
tions, has clearly shown the sensitivity of the results to both parts of the input of the calcula- 
tions. This is true not only in the Ip-lh responses, where the result was expected since they 
are dominated by the uncorrelated terms, but also in the 2p-2h responses, where the uncorre- 
lated term is not present. Physically meaningful results can be obtained only if single particle 
wave functions and short range correlations are linked through the nuclear hamiltonian. 

In our calculations the presence of short-range correlations does not produce remarkable 
effects. There are other effects beyond the pure mean-field model, such as final state inter- 
actions and meson exchange currents |^], which are more important. We should, however, 
point out that we have used only scalar correlation functions. There are indications that 
state dependent correlations together with MEC [^] could produce sizable effects. 

The next step of our work will be to apply the model to the description of exclusive 
electron scattering experiments, (e,e'N) and (e,e'2N), where one expects that the effects of 
the short-range correlations could be disentangled. 
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Appendix A. One-particle emission matrix elements 



In this appendix we give the expressions of the matrix elements involved in the one-nucleon 
emission part of the longitudinal and transverse responses. As discusses in Sect. ^, only 
the charge operator contributes to the longitudinal response, while both the electric and 
magnetic parts of the magnetization current, Eqs. ( p6| , p^ ) , contribute to transverse response. 
The convection current is considered in the one-point diagram only. 

We begin with the Coulomb (charge) matrix elements for which we have: 

^(ipih) ^ Gue)Kt.ipH^T\^)^ (76) 
Mir'^ = Gue)Ktu^-iiH^t^'^'\<i)^ (77) 

^S"'"^ = G'^(9')'^t.A7AE'^*..t.^(^^ + ^^ + ^)i' (78) 

iL 

( ji jh L y M[ih]. , 

Mir'^ = K,t,^n^Uj-I^ff'^^^'^\q), (79) 

^i.'f'^ = GW)K,t,7'pHT.^U,tr,(ilp + h+L)3f (80) 

iL 

( jp ji L V ^[ih];[pi], . 



ikL 

Jp 3k L \ ^[ih]f^^-^ ,j[pk];[ki] 
2 2 " 



^3.\"''^ = Gl{q^) St^,t, l^h E ^UJ, K,t, Sj^,^ + k+L) + k)3l (81) 

ikL 

Mi.r'^ = G|(,^)<^,,,,47r7,^.E^*..*.'^....e(^P + ^0xf (g) (82) 

i 

M'it''^ = G%{q^) St^,t, lih E ^ta, K,t, m + lh + L) + Ih) 11 (83) 



ikL 
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2 

jk jh L \ jH(^) j[kh];[ik] 



1 _i 

2 2 



Mi^s'''"^ = Ktu 4- E ^U,t,. K,t,. i-ir^^^^'^^ 7.i 7p^. (85) 

M^^t"-^ = K,t, 4vr ^ 7/. E G^i.('?') (7.i)^xf '(9) ^^f^''^^ ■ (86) 

In the previous equations we have indicated with G%{q'^) the electric form factor of the 
nucleon in the single particle state a, and we have defined the integrals 

lf\q) = j drr^jj{qr)Rl{r)Rp{r), (87) 



^ff'^'^k?) = j dnrljj[qr)Rl{n)Rp[n)-H^\n), 
X[pfh^\^q) = j dr^rljj{qr)p,^{n)nf\n), (90) 

where jj{x) is a spherical Bessel function, '^^^'^^(r) is given by Eq. (|4^), '^[(''(r) is given by 
Eq. (II), and 

/>ch(r) = Gr*(g') Pprot(r) + Gr'(g') Pneut(r) , (91) 

with Pprot('^) and Pncut('^) the proton and neutron densities, respectively. Finally, the integrals 
J are defined in Eqs. (Pg). 

The matrix elements of the transverse electric part of the magnetization current are given 

by: 

^E(lplh) _ G\i{q^) J Xp-Xh ^[pfe]/N I'qr.N 

^E(lplh) _ G'M(g^) , A^J Xp-Xh j[ph];[p]( . 
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P iL 



Ji Jh 



1 



L 




^=2 



K 



L K J 

Jp Jh Ji 



X 

1 K A 
1-10 



Ip li A 




K J L 
1-10 



1 



A 1 K 



(94) 



0, (95) 
9 %^ <^t,A (-l)'^+^-^+^'' J E '5*.,*. (-l)''^"^ C(Zp + k + L) 

P iL 

lijf ( I i\ ) Xg*'](g) Y.aL + J + \)X'(^ I; 0^ ) (96) 

j^-^ j L K J \ ( I K \\(k J l\\ ^; \ ^; 



c J Xp ~ Xh \ ^ AAA 



a/p + /. + i^)a/p + /.)?! (I \ yxfl(g)jif^]^['=^ 
-9 ^ K.th 7^/1 J E ^ti.tp <^Ji.> '^(^P + 

2M '''^'^ -yj J[J + 1) •^'^ 



jfeL 

2 



(97) 



(98) 



(99) 



2 2 
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^E(lplh) 
-'3.4;m 



(100) 



^E(lplh) 
3.5;m 



2M„ 



E (-i)^-^+^'=+^+^ (101) 



^E(lplh) 
3.6;m 



J \2 
ik) 



(102) 



In these equations we have indicated with G^(g^) the magnetic form factor of the nucleon 
in the single particle state a, and with Mq. The symbol Xa has been defined as: 



Xa = {la - ja) {"^ja + 1) • 



(103) 



For the transverse magnetic part of the magnetization current the corresponding matrix 
elements are: 



yM(ipih) _ _. Gljifl J 1 



s=±l 



(104) 



^M(lplh) _ _ 



(105) 



^M(lplh) 
-'2.2:m 



(106) 
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^M(lplh) 
-'2.3:m 



^M(lplh) 
-'2.4;m 



^M(lplh) 
-'3.1;m 



^M(lplh) 
-'3.2:m 



^M(lplh) 



L J' X 
1 K J 



s=±l 

Jp Jh Jt j 

0, 

2Mp V 27r ^ ^ 

iL 



2 \ Ip h 




L J' X 




2 ^'P 

I- i i- 
A 1 



(107) 



s=±i 



1 J 





L 


1 

2 














f k 











(108) 



L J' X 




A 1 K 



"^^ "oXf ^Ph- 2^ 



c(Zp+^.+i^)e(^p+yii( t !_\ ) jr'^t'^^ 



(109) 



E * \/ "^j^^s ^ + X/i + s J + 5s,i) (g) , 



s=±l 



2M„ 



■ (^tp,t^ 47r rj^h -J- E ^ti,tp <^jijp C{lp + ^i) ^0 



H;[p] 



(110) 



ikL 
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s=±l 



(111) 



^M(lplh) 
-'3.4;m 



-%q 



[ih];[p] 



s=±l 



(112) 



^M(lplh) 
3.5;m 



2M„ 



-^tpA 4^ 1? E ^t.-t. (-l)^->+^-'=+^+i 4^ 



J2 f 



1 
L2 



pi\;[kh] 



L 



(113) 



^M(lplh) 
-'3.6;m 



0. 



(114) 



In the previous equations we have defined 



-^(-l)^''+Me(/„ + /^ + A+l)iJ;jA (^-^f \ 



and J' stands for J + s. 

Finally, for the transverse electric and magnetic parts of the convection current contribut- 
ing to the one-point diagram we have: 

^E(lplh) ^ 1 gg(g^) c J 1 

q 2Mp 

{[(Xp - X^)(Xp + X^ + 1) + J{J + m^'^'^Q) + (115) 
[{Xp - Xh){Xp + Xh + 1)- J{J + 1)] Tj"'""^ {q) 



^M(lplh) 



2M„ 



V^J + l) 
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where we have used the fohowing definitions: 

= ldrrjj{qr)R;{r)^^, (117) 

= Jdrrjj{qr)^^Rh{r), (118) 

lfP''\q) = j dTTjj{qT)Rl{T)Rn{r). (119) 
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Appendix B. Two-particle emission matrix elements 

In this appendix we present the expressions of the matrix elements involved in the two-nucleon 
emission part of the longitudinal and transverse responses. The charge operator contributes 
to the first one and the corresponding Coulomb matrix elements are: 



(2p2h) _ 



2.1 



g'e (q') K Ai K a. ^ (-i)"" JpJhJT.T ^v-^^ ^p-f^- 

LK ^ 



jpi jp2 '^P 

jhi jh2 '^h 

K L J 



]-^JL ' 



(120) 



iL ^ 



J j jp2 jh2 L [ j Jp Jh J 



(121) 



jhi ji "^h \ \ ji jpi jp2 



^(2p2h) ^ Gf(g2)5,,^,,^5,,^,,^4;r(-l)^+^-+^''.+^J,J,5:5,,,,^ 



iL 



(122) 



.,J I jhi jpi L [ \ Jh Jp J ] ^[ih2], ^ ^\pihi];\p2i] 



Ti/i2 7p2i ^Pihi 



jp2 ji "^P 



Ji Jh2 Jhi 



M. 



(2p2h) 



3.9 



h2 



(123) 



Jhi 



Jh J "^P { ■j-[pih2] 
Jpi JP2 Jh2 



I if^'^^\q) j2 5,,,^^ {-ir^^^. ^,iki7^2iJt'^''^'''^ 



M. 



(2p2h) 



3.11 



G''E{<f)K^t,^ K2^t.l ^^.2JHi ^TT {-iy+^^l+^^2+^ aip2 + Ihl) JpJh (124) 



'^Pih2 



Jh J Jp 
Jpi Jhi Jh2 



jbi'iz]^^^ j^2hi];[p\ 



As already stated, in our model, the convection current is considered only in the uncorre- 
lated diagrams which do not appear in this case. As a consequence, only the magnetization 
current contributes to the electric and magnetic transverse terms. For the electric ones, the 
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corresponding matrix elements can be written as follows: 



^E(2p2h) 
-'2.1;m 



2Mp, 



K 



^( K J L 
1 -1 



J Pi JP2 

K L 

I 1 

'pi 2 

A 1 



^E(2p2h) 
-'3.1;m 



2M, 



pi 



■Ki,tn, K..tK, V^i-y+'^^+^'^^+'^^+'jpJhJjpi 



E mpi +k+j)Y-, 



r\p2h2\;[ihi\ 



iL 



Jp Jh J 
Ji Jpi Jp2 



L2 

jp2 i/t2 

jhi ji Jh 



ji jpi J 

_i _1 1 
2 2 



^E(2p2h) 
3.5;m 



2M 



4Tr{-iy+^^^+^^2+^JpJhJjh, 



^P2 



T\p2i\;\p\h\\ 



iL 



J h Jp J 
Ji Jh2 Jhi 



L2 

jhi jpi L 
jp2 ji Jp 



jh2 ji J 



1 1 



^E(2p2h) 
-'3.9;m 



2M, 



pi 



K2,tu, ^2,^., V4^ i-iy+^^i^^-i^^^ 



^.{Ipi + lh2 + J) Jp Jh J jpi 



jh2 \ Jh J Jp 

iZ.^ \ -^Pi ■^'»2 



jh2 jpi J 
2 2 



E ^u,t.2 (-i)^'-""^'"^ ^7.i7p^,. Ji^^'^f^'^' 



^E(2p2h) 
-'a.llim 



2M, 



pi 
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(.{Ipi + lh2 + J) C{lp2 + ^/ll) Jp Jh Jjpi jh2 

{Jh J '^P [ I 3h2 jpi J \ 7-b'l'*2] rr\p2h\];[p] 
JV. JH2 -\ 1 J ^ ^'^'^^ 

Finally, the magnetic transverse matrix elements of the magnetization current are: 



(130) 



Pi 



LK ^ 



^pi 2 ^Pi 

lh\ 2 jhi 

A 1 K 




jpi jp2 Jp 

jhi J/12 "^h 

K L J 



'Pi ^ ^hi 





L J' X 
I K J 



L J' X 




2Mp, 

jp2 3h2 ^ 
jhi ji 'hi 



J„ J) 



A. A A I nhi Ip2h2 L 
Jl Jpi Jp2 I 



p '^h 



E 



U + 5, 



s,-l 



J + 5^ 



•s,l 



(131) 



pM(2p2h) 
3.5;m 



= 4- (-i)-.+>.+^ i 



jhi jpi L 
jp2 ji Jp 



J2 



iL 



L2 



Jh Jn 



E 



8,1 



ife jh, J ^fi'^i^f^^^-^ 
-(x/.2 + Xi + ^^ + M4'''k«)' 



(132) 



^M(2p2h) 
-'3.9;m 



2M 



pi 



?^/i2 ^P2 '^^ 



Jpi Jp2 J/l2 



E 



s,l 



-{Xpi+Xh2+sJ + S,,i)I^P'''\q) 



(133) 
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\p2i];lihi] 



iL 



■,M(2p2h) 
3.11;m 








[P2/ii];[p] 





(134) 



The integrals and symbols used in these equations are defined in Appendix A. 
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